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Abstract 

We give a brief review of recent QCD analysis [I] carried out over the 
deep inelastic scattering data on F2 structure function and in the non- 
singlet approximation to the accuracy up to next-to-next-to-leading- 
order. Specifically, analysis was performed over high statistics deep 
inelastic scattering data provided by BCDMS, SLAC, NMC and BFP 
collaborations. For the coupling constant the following value a s (M|) = 
0.1167 ± 0.0022 was found. 



1 Introduction 

It goes without saying how it is crucial to know as accurate as possible 
the parton distribution functions (PDFs) and the value of the strong 
coupling constant in order to be able to make (relatively) solid predic- 
tions for various processes studied in a number of experiments. Within 
this realm, the deep inelastic scattering (DIS) of leptons off hadrons 
serves to be a cornerstone process to study PDFs which are universal 
and feed them further to other processes. 

Nowadays the accuracy of data for DIS structure functions (SFs) 
makes it possible to study Q 2 -dependence of logarithmic QCD-inspired 
corrections and those of power-like (non-perturbative) nature in a sepa- 
rate way (see for instance [2] and references therein) which is important 
for the analysis to be performed according to a well defined scheme. 
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Until recently a commonly adopted benchmark tool for the analysis 
happened to be there at the next-to-leading-order (NLO) level. However 
there have already appeared papers in which QCD analysis of DIS SFs 
has been carried out up to the next-to-next-to-leading order (NNLO) 
(see pQ and [3]-|10| and references therein). 

In the present paper we show the results of our recent analysis pQ of 
of DIS SF F 2 (x, Q 2 ) with SLAC, NMC, BCDMS and BFP experimental 
data involved [H]-[14] at NNLO of massless perturbative QCD. 

As in our previous paper the function F2(x,Q 2 ) is represented as a 
sum of the leading twist F2® CD (x,Q 2 ) and the twist four terms: 

F,( X ,Q 2 ) = F^ GD ( X ,Q 2 )(l + ~^\ . (1) 

As is known there are at least two ways to perform QCD analysis over 
DIS data: the first one (see e.g. [151 05]) deals with Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) integro-differential equations |17j and 
let the data be examined directly, whereas the second one involves the 
SF moments and permits performing an analysis in analytic form as 
opposed to the former option. In this work we take on the way in- 
between these two latter, i.e. analysis is carried out over the moments 
of SF F| (x,Q 2 ) defined as follows 

MpQCD/tviisty...^ = f 1 x n-2 p pQC D / t W ist2 / . . . ^ q2) ^ (2) 

JO 

and then reconstruct SF for each Q 2 by using Jacobi polynomial expan- 
sion method [H]-[20j (for further details see [21] ). The theoretical input 
can be found in the papers |21[ I22j. 

2 A fitting procedure 

To cut short this follows along the lines described in the paper |21j . 
Here we just recall some aspects of the so-called polynomial expansion 
method. The latter was first proposed in [23] and further developed 
in [23]. In these papers the method was based on the Bernstein polyno- 
mials and subsequently used to analyze data at NLO [25j[26] and NNLO 
level [HE]. The Jacobi polynomials for that purpose were first proposed 
and then subsequently developed in [181 EH 120] and used in [3]- [6], [27] . 

With the QCD expressions for the Mellin moments M^{Q 2 ) analyti- 
cally calculated according to the formulae given above the SFF 2 fc (x,Q 2 ) 
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is reconstructed by using the Jacobi polynomial expansion method: 



F*(x,Q 2 



x a (l 



Nmax 

' E 

n=0 



a,P)M* +2 (Q 2 



j=0 



where 0^' 6 are the Jacobi polynomials and a, 6 are the parameters fitted. 
A condition put on the latter is the requirement of the error minimiza- 
tion while reconstructing the structure functions. 

Since a twist expansion starts to be applicable only above Q 2 ~ 1 
GeV 2 the cut Q 2 > 1 GeV 2 on data is applied throughout. 

MINUIT program [28j is used to minimize two variables 



XSF 



pexp 



pteor 
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Aslope 
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where D = dlni^/dlnlnQ 2 . Quality of the fits is characterized by 
X 2 /DOF for the SF F 2 . Analysis is also performed for the SF slope D 
that serves the purpose of checking the properties of fits. 

We use free normalizations of the data for different experiments. For 
a reference set, the most stable deuterium BCDMS data at the value of 
the beam initial energy Eq = 200 GeV is used. 



3 Results 

In the nonsinglet approximation gluons are not taken into account, 
hence the cut on Bjorken variable (x > 0.25) imposed where the gluon 
density is believed to be negligible. 

The starting point of the evolution is taken to be Qq = 90 GeV 2 . 
This latter value is close to the average value of Q 2 spanning the cor- 
responding data. On grounds of previous knowledge the maximal value 
of the number of moments to be accounted for is N max = 8 [HI [20] 
(though we check N max dependence just like in the NLO analysis) and 
the cut 0.25 < x < 0.8 is imposed everywhere. 

Analogously to what has been done in [211 [T] the cut on the Bjorken 
variable x is imposed in combination with those placed on the y variable 
as follows: 



y 


> 0.14 


for 


0.3 < x < 0.4 


V 


> 0.16 


for 


0.4 < x < 0.5 


y 


> 0.23 


for 


0.5 < x < 0.6 


y 


> 0.24 


for 


0.6 < x < 0.7 


y 


> 0.25 


for 


0.7 < x < 0.8 
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which are meant to cut out those points suffering large systematic errors. 
Then a full set of data to be analyzed consists of 797 points. 

To verify a range of applicability of perturbative QCD we start with 
analyzing the data without a contribution coming from the twist-four 
terms (which means F2 = F^® 00 ) and perform several fits with the cut 
Q 2 > Qmin gradually increased. From Table 1 it is seen that unlike 
the NLO analysis quality of the fits starts to appear fairly good from 
Q 2 = 7 GeV 2 onwards. For the purpose of comparison, the twist-four 
corrections are added and the data with a general cut Q 2 > 1 GeV 2 
imposed upon is fitted. It is clearly seen that as in the NLO case 
here the higher twists do sizably improve the quality of the fit, with 
insignificant discrepancy in the values of the coupling constant to be 
quoted below. 

Table 1. a s (M§) and \ 2 i n the case of the combined analysis 



o 2 


N of 
points 


HTC 


X 2 (F 2 )/DOF 


a s (90 GeV 2 ) ± stat 


a s (M 2 ) 


1.0 


797 


No 


2.20 


0.1767 ± 0.0008 


0.1164 


2.0 


772 


No 


1.14 


0.1760 ± 0.0007 


0.1162 


3.0 


745 


No 


0.97 


0.1788 ± 0.0008 


0.1173 


4.0 


723 


No 


0.92 


0.1789 ± 0.0009 


0.1174 


5.0 


703 


No 


0.92 


0.1793 ± 0.0010 


0.1176 


6.0 


677 


No 


0.92 


0.1793 ± 0.0012 


0.1176 


7.0 


650 


No 


0.92 


0.1782 ± 0.0015 


0.1171 


8.0 


632 


No 


0.93 


0.1773 ± 0.0018 


0.1167 


9.0 


613 


No 


0.93 


0.1764 ± 0.0022 


0.1163 


10.0 


602 


No 


0.92 


0.1742 ± 0.0023 


0.1154 


11.0 


588 


No 


0.91 


0.1718 ± 0.0027 


0.1144 


12.0 


574 


No 


0.92 


0.1717 ± 0.0029 


0.1143 


13.0 


570 


No 


0.92 


0.1710 ± 0.0030 


0.1140 


14.0 


562 


No 


0.92 


0.1712 ± 0.0032 


0.1141 


15.0 


550 


No 


0.91 


0.1715 ± 0.0033 


0.1142 


1.0 


797 


Yes 


0.98 


0.1772 ± 0.0027 


0.1167 



The following values for parameters in the parametrizations of the 
parton distributions for the case corresponding to the last row of the 
above table are obtained: 

A^ s = 2.54, A% s = 2.38, A C NS = 3.29, A F N e s = 2.35, 
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b H N % = 4.16, b^ s = 4.22, b^s = 4.23, b F N e s = 4.39, 
d% a s = 6.08, d° 2 s = 3.89, d% s = 2.02, d^ e s = 3.31. 

The parameter values of the twist-four term are presented in Table 2. 
Note that these for H2 and D2 targets are obtained in separate fits 
by analyzing SLAC, NMC and BCDMS datasets taken together. It is 
seen that the values at NLO and NNLO match within errors with an 
average value being slightly less for the latter. Note that the values of 
this parameter at NLO is still compatible with those of the analogous 
analysis carried out in |15| . 

Table 2. Parameter values of the twist-four term in NNLO analysis 



X 


Ha (x) of H2 ± stat 


h^{x) of D2 ± stat 


0.275 


-0.183 ± 0.020 


-0.197 ± 0.009 


0.35 


-0.149 ± 0.028 


-0.171 ± 0.015 


0.45 


-0.182 ± 0.029 


-0.033 ± 0.031 


0.55 


-0.236 ± 0.052 


0.142 ± 0.057 


0.65 


-0.180 ± 0.135 


0.295 ± 0.108 


0.75 


-0.177 ± 0.182 


0.303 ± 0.158 



For illustrative purposes we visualize them in Fig. 1, where fairly 
good agreement between higher twist corrections obtained at NLO and 
NNLO is observed, that is not in contradiction with earlier studies (see, 
for example, [8]). 

We would like to note that the cut of the BCDMS data, which has 
increased the a s values (see Fig. 2 in pQ ) essentially improves agreement 
between perturbative QCD and experimental data. 

Indeed, the HTCs that are nothing else but the difference between 
the twist-two approximation (i.e. pure perturbative QCD contribution) 
and the experimental data are seen to become considerably smaller at 
NLO and NNLO levels, as compared to the NLO HT terms obtained 
in [15] and also with the results of analysis obtained within a fixed- 
flavor-number scheme (with a number of flavors fixed to be 4 and no 
F-cuts imposed on the BCDMS data (see Fig. 2)). 

Thus, using the analyses based on the nonsinglet evolution of the 
SLAC, NMC, BCDMS and BFP experimental data for SF F 2 with no 
account for the twist-four corrections and the cut Q 2 > 8 GeV 2 imposed, 
we obtain (with x 2 /DOF = 0.93) 

a s (Mf) = 0.1167 ± 0.0008 (stat) ± 0.0018 (syst) ± 0.0007 (norm) (3) 
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Figure 1: Comparison of the HTC parameter h±(x) obtained at LO, NLO and 
NNLO in the case of hydrogen data (the bars denote statistical errors). 



or 

a s {M 2 z ) = 0.1167 ±0.0021 (total exp.error). (4) 

Upon including the twist-four corrections, and imposing the cut 
Q 2 > 1 GeV 2 , the following result is found (with x 2 /DOF = 0.98): 

a a (Mz) = 0.1167 ± 0.0010 (stat) ±0.0020 (syst) ±0.0005 (norm) (5) 

or 

a s [M 2 z) = 0.1167 ±0.0022 (total exp.error) (6) 



4 Conclusions 

In this work the Jacobi polynomial expansion method developed in [181 
[19j [20] was used to perform analysis of Q 2 -evolution of DIS structure 
function F% by fitting respective fixed-target experimental data that 
satisfy the cut x > 0.25. Based on the results of fitting the value of 
the QCD coupling constant at the normalization point was evaluated. 
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Figure 2: Comparison of the HTC parameter h^x) obtained at LO, NLO 
and NNLO in the case of hydrogen data within a fixed-fJavor-number scheme 
(rif — 4) and with no y cuts imposed on the BCDMS data. 



Starting with the reanalysis of BCDMS data by cutting off points with 
large systematic errors it was shown pQ that the values of a s (M|) rise 
sharply with the cuts on systematics imposed. The values a s (M§) ob- 
tained in various fits are in agreement with each other. An outcome is 
that quite a similar result for a s (M z ) was obtained PQ in the analysis 
performed over BCDMS data (with the cuts on systematics) and that 
derived in the analyses done over the data of the rest, thus permitting 
us to fit available data altogether. 

It turns out that for Q 2 > 3 GeV 2 the formulae of pure perturbative 
QCD (i.e. twist-two approximation along with the target mass correc- 
tions) are enough to achieve good agreement with all the data analyzed. 
The reference result is then found to be 

a s {M 2 z ) = 0.1167 ± 0.0008 (stat) ± 0.0018 (syst) ± 0.0007 (norm), (7) 

Upon adding twist-four corrections, fairly good agreement between 
QCD (i.e. first two coefficients of Wilson expansion) and the data start- 
ing already at Q 2 = 1 GeV 2 , where the Wilson expansion begins to be 



7 



applicable, is observed. This way we obtain for the coupling constant 
at Z mass peak: 

a s {M 2 z ) = 0.1167 ± 0.0007 (stat) ± 0.0020 (syst) ± 0.0005 (norm) . (8) 

Note that there also is good agreement with the analysis |29| of the 
combined HI and BCDMS data, which was published by HI collabora- 
tion. Our result for a s (M z ) is also compatible with the world average 
value for the coupling constant presented in the review [30] , 

a s {M 2 z ) = 0.1184 ± 0.0007. (9) 

Concerning the contributions of higher twist corrections in the present 
work the well-known x-shape of the twist-four corrections while going 
from intermediate to large values of the Bjorken variable x is well re- 
produced. 

As a next step, we plan to study the DIS experimental data in the 
non-singlet approximation along with an application of resummation- 
like Grunberg effective charge method [31 j (as it was done in [27] at the 
NLO approximation) and the "frozen" |32L 133] and analytic [34] versions 
of the strong coupling constant (see [33j [35j [36] for recent studies in this 
direction). While our less immediate goal, a singlet analysis combined 
with a non-singlet one, is intended to be considered, with the effects 
of heavy-quark thresholds taken into account, as a follow-up to the 
previous study [37] that deals with nonsinglet part only. 

The work was supported by RFBR grant No.l0-02-01259-a. The 
work of GP was supported by the grant Ministerio de Ciencia e Inovacion 
FPA2008-01177. 
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